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9.1 Introduction 

Techniques previously designed for electroencephalographic (EEG) seizure detection 
in the newborn have been relatively inefficient due to their incorrect assumption 
of local stationarity of the EEG. To overcome the problem raised by the proven 
nonstationarity of the EEG signal, current methods are extended to a time-frequency 
(TF) approach [8, 10]. This allows the analysis and characterization of the different 
newborn EEG patterns, the first step toward an automatic TF seizure detection and 
classification. An in-depth analysis of the previously proposed autocorrelation and 
spectrum seizure detection techniques identified the detection criteria that can be 
readily extended to the TF domain. We present the various patterns of observed 
TF seizure signals and relate them to current specialist knowledge of seizures. In 
particular, initial results indicate that a quasilinear instantaneous frequency (IF) can 
be used as a critical feature of the EEG seizure characteristics. These findings led to 
propose a TF-based seizure detector. This detector performs a two-dimensional (2D) 
correction between the EEG signal and a reference template selected as a model of 
the EEG seizure in TF domain. 
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9.1.1 Seizures and the electroencephalographic signal 

Seizure may be the most frequent, and often the only, clinical sign of central nervous 
system dysfunction in the neonate [34], Unlike adult seizure, the clinical manifest- 
ations of newborn seizure are subtle and hence require the constant attention of a 
medical specialist for diagnosis. 

EEG has become a successful means of seizure detection in adults. This usually 
involves identifying sharp repetitive waveforms that indicate the onset of seizure. In 
adults, these EEG seizures are easily recognizable against a low amplitude random 
background characteristic of normal brain activity. The problem of detecting seizure 
in newborn babies, however, is complicated by a number of factors [28]. First, 
healthy newborn EEG signals representing normal brain activity often contain patterns 
such as spurious waveforms and sharp spikes. These characteristics, which would 
otherwise be detected as seizure in adults, are simply the result of extra electrical 
activity produced by the immature brain as it continues to form. Seizures, however, 
still usually appear in the EEG data as repetitive waveforms and the problem lies 
in discerning the healthy spikes from those formed from seizures. Second, visual 
symptoms of seizure, such as muscle spasms, rapid eye movement and drooling, are 
much subtler in newborns and may be easily missed. These visual indicators are 
also natural movements common to all newborn babies. Third, physical activity of 
babies in the intensive care environment is often subdued by medication to prevent 
injuries caused by unpredictable movements. This also reduces the chance of seizure 
detection using visual signs altogether. For these reasons, the EEG may be the only 
alternative for seizure detection in newborn. 



9.1.2 Nonstationary and multicomponent characteristics of the 
electroencephalogram 

There are three published methods for EEG seizure detection in newborns. The 
technique of Roessgen, Zoubir and Boashash [28] is a parametric approach based on 
a nonlinear estimation of 1 1 model parameters for detection. The two other methods 
are nonparametric. The spectrum technique introduced by Gotman et al. [16] uses 
frequency analysis to determine the changes in the dominant peak of the spectrum 
of short epochs of EEG data. The autocorrelation technique introduced by Liu [20] 
performs analysis in the time domain using short epochs of EEG data. 

All three techniques are based on the assumption that the EEG signals are locally 
stationary. However, a closer examination of these signals clearly shows that they 
exhibit significant nonstationary and multicomponent features [9, 24] (Figure 9.1). 
Therefore, the basic assumption on which these three methods are based is essentially 
invalid and at best only an approximation. 

To take this fact into account, we extend the autocorrelation and spectrum tech- 
niques to the nonstationary case by introducing both time variation and frequency 
variation parameters in the basic equations on which these techniques are based. 
This extension leads to a TF formulation in the (t - f) domain or a time- varying 
autocorrelation formulation in the time-lag domain (t — r). 




Time-Frequency Methodology for Newborn EEG Seizure Detection 



341 



Fs=20Hz N=600 



Time-res=5 




Frequency (Hz) 




FIGURE 9.1 

Time-frequency (TF) representations of newborn EEG seizure signal using the B distribution. This EEG 
clearly shows the multicomponent and nonstationary nature of the seizure. 



Having taken into account the nonstationary aspect of the EEG signals, we still need 
to deal with the other aspect, namely, the multicomponent behavior. This requires 
a selection of an appropriate time-frequency distribution (TFD) that is capable of 
handling multicomponent signals. Once this has been done, a calibration process 
is undertaken to translate the seizure criteria from the spectrum and autocorrelation 
domains to the TF domain. This initially involves the reproduction of the seizure 
detection criteria used by the two previously mentioned methods and map them in a 
joint TF domain. Features in the joint t- f domain characterizing seizure are then 
identified and a detection process based on those features is constructed and tested. 
Once we have a better understanding of the processes involved for each method and 
determined which components are most suited to a TF representation, we combine 
them in a TF-based detection process as illustrated in Figure 9.2. 

9.1.3 Electroencephalographic acquisition 

The electrical signals produced in the brain are monitored in a noninvasive manner 
by measuring variations in potential on the scalp. The EEG signals are measured 
using strategically placed small electrodes on the scalp. One electrode, usually at the 
base of the skull, acts as a reference (ground) signal, and various channels of data are 
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FIGURE 9.2 

Tt-bascd seizure detection process. 



created by measuring the voltage dilfercnces between neighboring electrodes. Due to 
the size of most newborn babies’ heads, only 5 channels of EEG have been recorded 
m each session using the 10 to 20 International System of Electrode Placement. The 
sampling rate used for recording the EEG data was 256 Hz. For the sake of artifact 
detection, three auxiliary signals representing electro-oculogram (EOG), electrocar- 
diogram (ECG) and respiration were also recorded. Data used in this study were 

collected at the Royal Women's Hospital Perinatal Intensive Care Unit in Brisbane 
Australia. 




9.2 Seizure Criteria and Time-Frequency Methodology 
9.2.1 Basic definitions and relationships 

The autocorrelation function of a random nonstationary process x(t) may be written 

nc* 



Rjc(E T) — k [x(f )j'(t + 7")] 



(9.1) 
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where r represents a time difference, or lag and E is the expectation operator. If the 
process is at least wide-sense stationary (WSS), the autocorrelation function becomes 
simply a function of r. that is: 

7?*(f,r)=f? r (r) (9.2) 

For the case of stationary signals. Wiener-Khintchine theorem states that the power 
spectral density of x(t) can be obtained by simply taking the Fourier transform (FT) 
of Equation (9.2). In other words: 

S r (/)= FT{R x (r)} (9.3) 

T-*f 

where R. c (r) is defined in Equation (9.2). This equation is valid only for stationary 
signals. For the case of nonstationary signals. Equation (9.3) gives averaged values 
over the variations of the signal spectral contents and thus significant information 
may be lost. 

For nonstationary signals, the Wiener-Khintchine theorem can be extended so that 
the time-varying power spectral density S x (t.f) of a signal x (I ) is also related to the 
time-varying autocorrelation function /? x (f.r) of the same signal by the following 
FT relation: 

S x (t.f)= FT {R x (t. r)} (9.4) 

T~*f 

or more explicitly: 

OC 

S r \tf)= I RAtr)e- j2nfT dT (9.5) 

— OO 



The expression 5 X ((. /) given by Equation (9.5) reduces to the Wigner-Ville spec- 
trum [5]. For practical reasons. r(t) is replaced by its analytical associate z(1). 
Estimates of S z (t,f) can be expressed by the quadratic class of TFDs [5 J: 

Pt(t,f) = W g (t,f)**7(t,f) (9.6) 

where ~>(t, /) is a 2-D kernel window that is application dependent. f ) is the 
Wigner-Ville distribution (WVD) of z(t ), and ** indicates a double convolution in 

both time and frequency. By expanding Equation (9.6). we obtain [5]: 



OC OO OG 



Pz(Ef) = III i/(Kr)z(u + r/2)f(u-r/2)( j2l(l ' 1 t/u * 'dududr (9.7) 



— OC — OC — OO 



where the superscript ’ stands for complex conjugate. The kernel filter q(o, r) char- 
acterizes a particular TF distribution and is related to y(f, /) via a double Fourier 
transformation. Equation (9.7) is the general formulation of the quadratic class of 
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FIGURE 9.3 

TF representation of a three-component signal. 



TFDs in terms of the Doppler-lag kernel g{v, r) [6]. The TFD p z (t, f ) can be further 
expressed in the following form: 

OO OO 

p z (t,f)= j J G(t — u,r)z(u + T/2)z*(u — T/2)e~ : ’ 27r f T dudT (9.8) 

— OO — OO 

where G(t, r) is the inverse FT of filter kernel g(y, r). A shorter version of Equation 
(9.8) is given by: 

p,(t,/)=FT{K,(t,T)'G(t,T)} (9.9) 

T~+f 

wher eK z (t,r) — z(t + r/2)z*(t — r/2) is the signal kernel. This equation leads to 
the easiest discrete-time version of TFDs for implementation purpose [5]. 

A typical TF plot of a nonstationary multicomponent frequency-modulated (FM) 
signal is shown in Figure 9.3. The representation on the left is the time-domain 
representation. From this representation we observe that the signal significantly 
changes its periodicity characteristics at time t = 100 sec and time t = 175 sec. 
The representation at the bottom is the magnitude spectrum of the signal. It shows 
that this signal consists of three disjoint frequency bands. The main plot is the TF 
representation of the same signal. It not only combines the information from both time 
and frequency representations but also preserves the phase information, disregarded 
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or unavailable in the other two representations. This allows the joint TF representation 
to clearly show that the signal consists of three components characterized by: 

1 . Frequency 0.05 to 0. 1 Hz during time t = 0 to 100 sec 

2. Frequency 0.25 to 0.3 Hz during time t = 100 to 175 sec 

3. Frequency 0.4 to 0.45 Hz during time t = 175 to 250 sec 

The TF representation shows not only the number of components present in the signal 
but also how they vary in time and frequency, when they start and finish and whether 
they coexist at the same time or same frequency. For this particular example, the TF 
representation indicates the linear FM nature of all three components, which cannot 
be deduced from either time or frequency representations only. 



9.2.2 Time-frequency distribution selection 

The EEG data collected thus far show that neonatal EEG seizures are highly nonsta- 
tionary and occasionally multicomponent, and are mostly concentrated in the band 
of frequency (0 to 5) Hz [8-10]. These factors must be taken into account when 
selecting an optimal TF distribution, because each TFD is more suited to representing 
signals with particular characteristics. For example, the WVD is known to be optimal 
for monocomponent linear FM signals but performs poorly for multicomponent oi 
highly nonlinear FM signals. 

Because neonatal EEG signals are nonstationary and occasionally multicompon- 
ent, the selected TFD needs to have good spectral resolution and reduced cross terms. 
For this sake, the performances of several distributions have been compared to find an 
optimal TF representation of real neonatal EEG data. The comparison study included 
the spectrogram, Wigner— Ville (WV), Choi— Williams (CW), B (BD), Zalto— Atlas- 
Marks (ZAM), Bom-Jordan (BJ) and Rihaczek-Margenau distributions [8]. The 
performances of the resulting TFDs have been compared using an objective quantita- 
tive measure criterion [7]. Based on this criterion along with extensive experimental 
results, the B distribution BD with the smoothing parameter (3 equals to 0.01 — 
see Equation (9.10) — has been found to be the most suitable TF representation of 
the newborn EEG. The BD is defined in terms of its time-lag kernel — see Equation 
(9.8) — which is given by: 



G(t,r) - 




( 9 . 10 ) 



where f3 (0 < (3 < 1) is a parameter that controls the sharpness of the cutoff of the 
2D filter in the time-lag domain ( t , r). Figure 9.4 shows the TFDs of a 30-sec 
sample of real newborn EEG data using the B, CW, ZAM and WV distributions and 
clearly illustrates the superiority of B distribution in terms of resolution and cross-term 
reduction. 
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ZAM distribution using parameter a- 2 WVD 

FIGURE 9.4 

TF representations of a real epoch of newborn EEG using different TFDs. 

9.2.3 From time domain to time-frequency domain 

Seizure detection using the autocorrelation method relies on the assumption that the 
essential characteristic in newborn seizure EEG is periodicity. The periodicity of the 
autocorrelation of short epochs of EEG data is quantified (scored) using a nonlinear 
rule-based algorithm. In this technique, an epoch consisting of 30 sec of data is 
divided into 5 windows (sec Figure 9.5). Depending on the autocorrelation function 
of each window, up to four primary periods (Tj, ... ,T 4) are calculated for each 
window in an epoch as shown in Figure 9.6. These times correspond to the moment 
centers of the lirst four peaks of the autocorrelation function. The windows are then 
scored whereby more evenly spaced primary periods are allocated larger scores. After 
each window in an epoch is scored, a second rule-based scheme is used to classify 
each epoch as positive or negative. If an epoch is classified as positive in two or more 
channels, this epoch is considered to contain seizure [20], 

9.2.3. 1 Seizure criteria for the autocorrelation method 

A single window is declared as seizure positive by the autocorrelation method if the 
following criteria are met [8]: 
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30 second epocr. 




FIGURE 9.5 

Epoch anil window definitions according to the autocorrelation method. 



1 At least four periods exist within the positive hall ol the autocorrelation 
function. 

2. The differences between the ratios of each moment center to the first and the 
nearest integer are less than 0. 1 50. 

3. The total score obtained by summing all moment center scores is greater than 
or equal to 12 (out of a maximum of 15). 



9.2.3.2 Mapping of autocorrelation seizure criteria to time-frequency domain 

An EEG signal within a given window is considered seizure if in the TF domain a 
continuous spectral peak exists within that window and meets the following criteria: 

1 . All frequencies within the spectral line are greater than 0.625 Hz within 
6.4-sec windows or greater than 0.909 Hz within 4.4-sec windows. 




13 14 15 16 17 18 19 



Time (s) 



FIGURE 9.6 

Autocorrelation function for one window. 
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2. The length of continuous dominant spectral peak within the window is greater 
than 3 sec. 

A continuous spectral peak is defined as adjoining peaks within the TF array above a 
selected threshold of one fifth the maximum array value. The threshold was chosen 
experimentally. 

Criterion 1 is a direct translation to the frequency domain of the first criterion 
described earlier in the autocorrelation domain. The second criterion has been 
deduced by observing several TF representations of seizure positive windows as 
defined by the autocorrelation method. The scoring system proposed by the auto- 
correlation method focuses on identifying periodic regions of data in the lag domain. 
Periodic regions are clearly identified in TF representations by a dominant spectral 
peak occurring for a certain time interval. 



9. 2.3.3 Results and discussion 

Extraction of the seizure criteria listed earlier in the TF domain has been successfully 
calibrated with the autocorrelation method. Peak detection techniques have been 
employed to simplify the extraction process, resulting in a detection array illustrating 
positions and lengths of continuous spectral lines within each epoch. Figure 9.7 
shows the algorithm flowchart used in the calibration of this method. 

Good results have been obtained using TF algorithm to detect individual seizure 
windows of real neonatal EEG in the TF domain. Approximately 75% of windows 
detected as seizure positive by the autocorrelation method are detected using TF 
criteria listed earlier. The results of applying the TF-based detection method are 
summarized in Figure 9.8. In this figure, the original epoch refers to the raw TF array 
produced from preprocessed EEG data. Image of this array appears on the left side 
of the figure. The image is also divided into four distinct 6.4-sec windows and one 
4.4-sec window as defined by the autocorrelation method. Window scores attributed 
to these windows are displayed at the end of each window division for easy comparison 

between the TF information and the corresponding score allocated by autocorrelation 
method. 



9.2.4 From frequency domain to time-frequency domain 

The spectrum analysis method has also been used to detect periodic discharges in 
the EEG signal. As shown in Figure 9.9, a background epoch is defined as a 20-sec 
segment of EEG finishing 60 sec before the start of the current 10-sec epoch under 
investigation. From the frequency spectrum of each 10-sec epoch, the following 
features are extracted: the frequency of the dominant spectral peak, the width of the 
dominant spectral peak and the ratio of the power in the dominant spectral peak to 
that of the background spectrum in the same frequency band. 



Time-Frequency Methodology for Newborn EEG Seizure Detection 



349 
















Applications in Time- Frequency Signal Processing 



Original Seizure Epoch 49 Part 1 Filtered Seizure Epoch 49 Part 1 

mm :iW, * £ 



v> -4 

«UU( I "o 

o 

u 

M 



8 

Mi n . 



-it .• I w 



12 3 4 

Frequency (Hz) 



0 

12 3 4 

Frequency (Hz) 



FIGURE 9.8 

The mapping of the features from autocorrelation lo TF. Seizure windows are defined by the presence of 
high value lines. As the length of each spectral line is denoted by its color, windows with dark lines should 
indicate scores higher or equal to 12. 
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FIGURE 9.9 

Epoch and w indow definitions according to the spectrum method. 
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*>.2.4.1 Seizure criteria in the spectral domain 

The 1 0-sec epoch of EEG data is considered seizure positive if any of the 3 following 
criteria are met [16]: 





Dominant 

Frequency 


Half- 

Maximum 

Bandwidth 


Power 

Ratio 


Criterion 1 


0.5- 1.5 Hz. 


<0.6 Hz 


3-4 


Criterion 2 


1.5-10 Hz 


<0.6 Hz. 


2-4 


Criterion 3 


1.5-10 Hz 


<1 Hz 


4-80 



To limit the number of false alarms, a seizure detection is discounted if the epoch is 
largely nonstat ionary, if a large amount of alternating current (AC) power noise is 
present or if it appears that an EEG lead has been disconnected. 

‘>.2.4.2 Mapping of spectrum criteria to time-frequency domain 

The preceding criteria pertaining to frequency and bandwidth are clearly observable 
in the TF domain. That is. each spectra containing a dominant peak that meets either 
of the criteria: 

1 . Frequency in the range 0.5 to 1 .5 Hz and width <0.6 Hz 

2. Frequency in the range 1.5 to 10 Hz and width <1 Hz 

may be considered for further seizure detection pertaining to power ratio. Disregard- 
ing the power ratio criteria, the second criterion given in the preceding table becomes 
a subset of the third criterion. 

*>.2.4.3 Results and discussion 

The calibration of the two TF criteria identified earlier is illustrated in Figure 9.10. 
This essentially extracts frequency and width information, the results of which are 
visible in the plots shown in Figure 9.1 1 for a 30-sec epoch. Data are presented by 
highlighting the dominant frequency with a color indicating the w idth of the spectral 
peak. Boxed sections of the array indicate EEG regions detected as containing seizure 
by the conventional spectrum method. 



9.2.5 Time-frequency patterns 

Calibration of the TF method using both autocorrelation and spectrum method has 
been successfully performed. This calibration led to consideration of a fully integrated 
TF detection method by combining the different TF seizure features. The visual 
analysis of these EEG TF features indicates that there are two distinguishable classes 
of TF EEG patterns: seizure and background. Further, the seizure patterns can 
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FIGURE 9.10 

Implementation and calibration of TF extension of frequency-based method. 



be characterized in terms of the dominant spectral peak as either a linear FM or a 
piecewise linear FM whereas the background patterns exhibit a low-frequency burst 
of activities or activities with no clearly defined pattern. These visual observations 
correlate well with the clinical information related to the different patterns found in 
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FIGURE 9.1 1 

Mapping of the features from spectrum to TF. The color of the dominant frequency is related to the width 
of the spectral peak. Boxed sections indicate the presence of seizure as delected by the spectrum method. 



EEG [2 1 1 . A few representative TF patterns that characterize each of the subclasses 
defined so far are given next. 

9.3 Seizure Patterns 

As mentioned earlier, the following patterns characterize the different patterns found 
thus far in the EEG signals in terms of dominant spectral peak. 

9.3.1 Linear frequency-modulated patterns 

A number of the EEG seizures analyzed in the TF domain can be approximated by 
linear FM (LFM) with either fixed or time- varying amplitudes. These patterns can 
be further classified into the follow ing subclasses: 

9.3. 1.1 Class A (linear frequency-modulated patterns with a quasi-constant 
frequency) 

Figure 9. 1 2 shows a seizure that has a linear FM behavior with an almost constant fre- 
quency. The amplitude of the TF seizure pattern increases at the onset and decreases 
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FIGURE 9.12 

TF representation of a seizure that shows a linear FM behavior with constant frequency. 



toward the end. These observations are in agreement with the seizure definition 
found in medical literature [21 ]. In this figure, we can clearly see the dominant peak 
(frequency domain) and the repetitive pattern of the spikes and sharp waves (time 
domain ) characterizing the newborn seizures. A major advantage of the TF represen- 
tation is that we can easily distinguish the seizure from the background. An interesting 
phenomenon that was observed is the existence of subharmonics in a number of TF 
representations of EHG seizure as illustrated in. for example. Figure 9.13. This 
characteristic of the newborn EEG seizure has also been reported by Lombroso [21 J. 

9.3. 1.2 Class It (short linear frequency-modulated patterns with a 
Quasi-constant frequency) 

Figure 9. 14 show s a seizure pattern that differs from those in class A in one important 
aspect, the duration. A disagreement exists between researchers about what consti- 
tutes a seizure. The duration of rhythmic discharges that characterizes a seizure is 
highly variable, from as short as 1 sec to as long as 30 min [25]. To consider an EEG 
signal as a seizure, some researchers require that it must last at least 10 sec 1 12. 16]: 
others require a minimum of 20 sec [23], whereas a third group does not specify a time 
limit 1 19-21 ]. Due to this disagreement, these controversial seizures were classified 
as a separate category. 
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FIGURE 9.13 

TF representation ot a seizure ihai shows subharmonics. 
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FIGURE, 9.14 

TF representation of a short duration harmonic discharge “seizure." 
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9.3.1.3 



Class C (linear frequency-modulated 
frequency) 



patterns with a decreasing 



Figure 9.15 of this class can be characterized by its IF that tends to decrease with 
time. This frequency-decreasing behavior is widely accepted as a characteristic of 
the EEG seizure [21, 29], The advantage of the TF representation over those of the 
time or frequency is that the linear FM character is easily recognizable. A similar 
behavior can be seen in the work of Franaszczuk on adult seizures originating from 
the mesial temporal lobe during the periods of organized rhythmic activity [14]. 



9.3.2 Piecewise linear frequency-modulated patterns 

A pattern of this class is shown in Figure 9.16. It can be approximated by piecewise 
linear FM. Something similar can be observed in one of the TF representations pro- 
duced by Franaszczuk in his study of the adult seizure originating from the mesial 
temporal lobe [14]. 
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FIGURE 9.15 

TF representation that shows an LFM behavior with decreasing frequency. 
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FIGURE 9.16 

TF representation of an EEG seizure that shows a piecewise LFM behavior. 



9.4 Electroencephalogram Background Patterns 

By background, we mean any signal that is not classified as seizure. Two distinct 
patterns have been observed while analyzing the EEG signals: burst suppression 
activity and irregular activity with no clear pattern. 

9.4.1 Class E (burst suppression) 

Figure 9.17 is an example of this class. It shows instance of burst activities. These are 
short period signals with a high energy lasting for few seconds and usually occurring 
at frequencies below 8 Hz. These features are characteristic of burst suppression as 
defined in the literature. Burst suppression consists of burst of high voltage activity 
lasting 1 to 10 sec and composed of various patterns (5 [0 to 4 Hz] and 6 [4 to 8 Hz] 
with superimposed and intermixed spikes, sharp waves and faster activity) followed 
by a marked background attenuation [29]. 

9.4.2 Class F (activities lacking a specific pattern) 

Figure 9. 1 8 is an example of EEG epochs lacking a well-defined and consistent pattern. 
These types of activities do not seem to be constrained within the low-frequency band 
characterizing the seizure signals. 
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FIGURE 9.17 

TF representation of burst suppression. 
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FIGURE 9.18 

TF representation of an EEG lacking specific and consistent pattern. 
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9.5 Time-Frequency Matched Detector 

The matched filter (which essentially reduces to a correlator receiver) is very popular 
for constructing detectors and classifiers. For known deterministic signals in white 
Gaussian noise (WGN), the correlator receiver is known to be optimal in the Neyman 
Pearson sense for detection [32] and in a Bayesian sense for classification [31]. To 
construct a detector, the output of a correlator is compared with a threshold. The 
threshold is chosen such that the probability of a false alarm is maintained. 

Traditionally, the correlator receiver is implemented in time domain as a one- 
dimensional (1-D) correlation between the received noisy signal x(t) and a reference 
signal s(t). For random signals, this can be done using the corresponding spectral 
representations. To extend this detector to handle nonstationary signals, the 1-D 
correlation is replaced by a 2-D correlation involving the TF distribution of x(t) and 
s(t), as shown in Figure 9.19. The resulting test statistic is given by: 

T(x) = J j p x (t,f)p s (t,f)dtdf (9.11) 

This type of detector has been implemented using different quadratic TF distributions 
such as the spectrogram [1, 33], the reassigned spectrogram [33], the WVD and cross- 
WVD [2-4, 13, 18, 30], and generalized time-shift covariant quadratic TFD [26-27]. 
It has also been implemented in the ambiguity domain by replacing the TFDs by the 
auto- and cross-ambiguity functions [33]. 

By using Moyal’s formula [17], we get: 

J j Px{t,f)p s {tJ)dtdf = J x(t)s*(t)dt (9.12) 

This equality is only valid when the filter kernel g(o, r) in Equation (9.8) is unimod- 
ular, that is: 

\g(u,T)\ = l (9.13) 

for all v and r. 




FIGURE 9.19 

TF matched receiver. 
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This is the case, for example, for the WV, the Rihaczek [13] and the generalized 
time-shift quadratic TFDs [27], Equation (9.12) is an alternative interpretation of the 
correlator receiver in terms of a correlation of the TFDs. 

It has been shown [2, 13] that for the case of deterministic signal in noise (even 
WGN) that the TF-based correlator is suboptimal. This suboptimality can be partly 
explained by the nonlinearity of the quadratic TFDs that accentuates the effects of 
noise by introducing artifacts. 

To use a correlator receiver, it is usually required that the wave shape of the reference 
signal (or other related information such as its TFD) as well as the noise statistics are 
known. For the case of EEG seizure detection, the wave shape of the EEG seizure 
is generally unknown. In Section 9.3, however, we have shown that the EEG seizure 
can be characterized by a linear or a piecewise linear FM. Based on this conclusion, 
we propose to construct a TF-based matched filter (i.e., essentially a correlator). To 
do this, a representative TF distribution of a linear FM p re f(t, /) is selected to serve 
as template (reference). The correlator statistic T(z ) used is the 2-D cross-correlation 
between the noisy signal and the reference signal TFDs [2], that is: 

T(z) = J j p re f(t,f)p* z (t,f)dtdf (9.14) 

where 2 is the analytical signal corresponding to the EEG signal to be analyzed.- 



9.5.1 Implementation of the time-frequency-matched detector 

A working model of the TF-matched detector as shown in Figure 9.20 has been 
implemented. A description of its main components are described in the next section. 

9.5.1. 1 Preprocessing 

This stage can include any preprocessing of the EEG signal collected from the new- 
born’s scalp. Some major reasons for preprocessing the EEG signal include artifact 
(such as ECG, EOG and EMG) removal, noise filtering and resampling the signal 
to comply with detector input specifications. In the present implementation of this 
detector, the process of artifact removal has not been included. 

9.5.1.2 Signal restructuring 

Because it is neither desirable nor practical to work with a very long EEG signal, we 
found it necessary to segment the EEG into an array of signals of fixed length. After 
some preliminary testing, we found a length of 2 min to be suitable for performing 
the cross-correlation. Shorter signal lengths led to higher rates of miss detections and 
false alarms. Once the full input EEG signal is divided into blocks of 2-min duration, 
each block is stored as a row of the newly formed array of signals. In the development 
of this detector, a protocol of 50% overlap of each block was adopted, as displayed 
in Figure 9.21. 
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FIGURE 9.20 

Block diagram of TF-matched detector. 




FIGURE 9.21 

Segmenting the EEG signal in blocks of 2-min length. 
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9.5.1.3 Detection loop 

The detection loop is executed until all the blocks of the input EEG signal have been 
processed. An offset value is maintained, giving a precise location in the original 
signal where abnormal events are detected. 

9.5.1.4 Cross-correlation implementation 

This forms the critical step in the effectiveness of this TF seizure detector. The cross- 
correlation between the input TF array and the template (mask) is obtained using the 
2-D cross-correlation function given by Equation (9.14). The most crucial process is 
the choice of the template. This problem is discussed later. We found it convenient 
to normalize TF distributions of both the reference signal and the EEG signal. 

9.5.1.5 Amplitude and length criteria 

Ideally there would be one peak value in the output of the cross-correlation array, 
with its amplitude determining the presence or absence of seizure. This proved a 
little unreliable, and it was decided to search for a sequential series of values over 
the amplitude threshold. This proved very successful, and a minimum ridge length 
of 20 sec over the amplitude threshold was classed as a seizure. The 20-sec length 
adopted is larger than the minimum 10-sec length of EEG seizure adopted by many 
neurologists [12]. This length can be changed at will. An example of the correlation 
output, following the application of a threshold is shown in Figure 9.22, where the 
ridges of length equal or greater than 20 sec indicate the presence of seizure. 




FIGURE 9.22 

Amplitude threshold applied to output of the matched filter. 
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FIGURE 9.23 

Different steps of seizure localization. 

9.5.1.6 Map seizure decision to real-time location 

This stage simply ties all the independent decisions on each block of processed signal 
(remapping any seizure decision to a time series function) of equivalent length to 
the input EEG signal. This output waveform consists of ones or zeros, where one 
indicates the presence of seizure at the corresponding time. Figure 9.23 summarizes 
the five stages from input signal to seizure indicated on the output signal. 

9.5.2 Experimental setup 




To validate and calibrate the detection method, simulated data generated by the EEG 
model proposed by Roessgen, Zoubir and Boashash [28] is used. 
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FIGURE 9.24 
EEG model. 



9.5.2.1 Electroencephalogram model 

Roessgen extended the model of the a rhythm of the thalamus proposed by da Silva 
[22] by introducing a seizure input waveform Z (t) as shown in Figure 9.24. The model 
considers a population of interconnected neurons driven by a random input — a WGN 
P(t) — that is assumed to originate from deeper brain structure such as thalamus and 
brain stem. The parameters of the model reflect some physiological characteristics 
such as neuronal interconnectivity, synaptic pulse response and excitation threshold. 
For an in-depth discussion of this model, see [28], 

To account for the nonstationarities of the EEG, the periodic sawtooth signal Z(t) 
of period f c used by Roessgen was replaced by the following more realistic quasi- 
periodic nonstationary signal that fits a linear FM law with a time- varying amplitude: 

S LFM (t) = Z(t)ei™ t2 (9.15) 

The typical range for a is [—0.07 to 0]. The initial frequency f 0 is typically in the 
range of 0.5 to 5 Hz. These values were reported in [15, 21]. A signal-to-background 
ratio (SBR) was introduced to account for the fact that the seizure and background 
components of the newborn EEG are of almost similar magnitudes. The SBR in 
decibels is defined as: 

SBR = lOlog (9.16) 

where a\ and represent the power in the seizure and the background components 
of the EEG signal, respectively. 
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FIGURE 9.25 
Template model. 



9.5.2.2 Template selection 

The reference template to be used for detection is the TF distribution of a linear FM 
(LFM, Figure 9.25). The duration of the LFM signal is set at 20 sec as discussed 
earlier. We have chosen the slope of the LFM to be negative because the majority of 
the seizures encountered at this stage are characterized by a frequency that is either 
constant or decreasing with time. 

To find the best slope of the template ridge (see Figure 9.25), that is, the one that 
corresponds to the best detection rate for the different types of seizures, a testing 
stage was necessary. A number of templates with different ridge slopes were created. 
These templates were tested using the seizure detection system presented earlier 
(see Figure 9.20) and a number of EEG signals produced by the model in Figure 9.24. 
These signals were generated using different combinations of the model parameters. 
It was found that the template with ( = —0.05 Hz/sec gives the best results. 

9.5.2.3 Threshold selection 

To select an optimum threshold that realizes the good compromise between the rate 
of good detections and the rate of false alarms, a similar test to the one previously 
discussed was performed. Among the different thresholds tested, the ones that have 
a value around 4 resulted in a good compromise. 
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9.5.3 Time-frequency detector evaluation 

To check its reliability and to validate its effectiveness, the TF detector is tested 
using synthetic data generated by the EEG model (see Figure 9.24). The synthetic 
signal contains a mixture of randomly placed seizures and background. For this 
investigation, 1000 synthetic EEG signals were used, with the parameter selection 
varied as follows: 

• SBR = random [0 to 20 dB] 

• SNR = random [0 to 20 dB] 

• a = random [—0.06 to 0] 

• f c = random [1.5 to —5] 

The average detection rate was 99. 1 whereas the false alarm rate was 0.4%. These 
good results were expected because the template adopted is well adapted to the EEG 
model that only produces noisy linear FM signals with time-varying amplitudes. 
Other EEG models could be able to generate different types of EEG-like seizures 
such as piecewise LFMs. For those signals, a template such as the one proposed by 
Celka, Boashash and Colditz [11] would be used. 



9.6 Discussion and Conclusions 

This chapter shows that the TF domain is the preferable basis from which to develop a 
complete EEG seizure detection scheme. The proposed TF approach was successfully 
calibrated using existing methods based on autocorrelation detection criteria and 
spectral estimation methods 

The initial patterns obtained by a TF analysis of EEG seizure signals specifically 
confirm that EEG seizures in newborn are well characterized by a linear FM or 
piecewise linear FMs. These results are encouraging and further analysis on other data 
sets is currently in progress to refine and extend these findings. The characterization 
of nonstationary EEG signals in the TF domain is the first step toward an automatic 
method of seizure detection and classification that use powerful tools of TF signal 
processing [6]. In this direction and based on the findings thus far obtained, we 
proposed a TF detector based on cross-correlating the TF representation of the EEG 
signal with a TF reference template. The design of the template takes into account the 
TF characteristics of the EEG seizure thus far extracted. The performance of this TF 
detector was tested on synthetic signals. The results obtained were very encouraging 
(99.1% detection rate and 0.4% false alarm rate). Current work not yet finalized at 
time of publication is concentrating on using real newborn EEG signals. 
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